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Abstract 

Degenerate pump-probe reflectivity experiments have been performed on a single crystal of 
bismuth telluride (Bi2Te3) as a function of sample temperature (3K to 296K) and pump intensity 
using ~ 50 femtosecond laser pulses with central photon energy of 1.57 eV. The time resolved 
reflectivity data show two coherently generated totally symmetric Ai s modes at 1.85 THz and 3.6 
THz at 296K which blue shift to 1.9 THz and 4.02 THz, respectively at 3K. At high photoexcited 
carrier density of ~ 1.7 x 10 21 cm -3 , the phonon mode at 4.02 THz is two orders of magnitude 
higher positively chirped (i.e the phonon time period decreases with increasing delay time between 
the pump and the probe pulses) than the lower frequency mode at 1.9 THz. The chirp parameter, 
f3 is shown to be inversely varying with temperature. The time evolution of these modes is studied 
using continuous wavelet transform of the time-resolved reflectivity data. 
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INTRODUCTION 



When a narrow band gap semiconductor is excited with intense femtosecond laser pulses, 
a dense electron-hole plasma is produced due to the promotion of electrons from bonding 
states to anti-bonding states, which can cause large changes in bond lengths leading to 
possible structural transitions . The photo-excitation of carriers changes the equilibrium 
positions of the atoms; the atoms then oscillate around their new equilibrium positions, a 
mechanism called as displacive excitation of coherent phonons (DECP) {s-jj. Thus, DECP 
is the dominant mechanism in opaque samples [3:] compared to impulsive stimulated Raman 
scattering (ISRS) j(| in transparent materials. Later, it was shown that DECP is a special 
case of ISRS when excited resonantly The investigations of coherent phonons performed 
under high photo-excited carrier density (PCD) > 10 20 cm -3 have been very few till now 



15] . In addi- 



In tellurium, time resolved reflectivity experiments [8] performed using 100 fs pulses at 
a PCD of ~ 5 x 10 21 cm -3 showed an instantaneous large red shift (13%) of Ai coherent 
phonon frequency, attributed to the electronic softening or bond weakening 
tion, the phonon time period decreasing with the delay time between the pump and probe 
pulses corresponds to an asymmetric line shape in frequency domain. This linear sweep 
in the frequency with the pump-probe time delay, termed as phonon chirping, originates 
from the rapid change of photoexcited carrier density across the sample thickness due to 



semimetal, bismuth 



carrier diffusion resulting in different amounts of phonon renormalization jl0|. In case of a 



12j,ll3| at carrier density of ~ 3 x 10 21 cm 3 in a degenerate pump-probe 



reflectivity experiment done at room temperature, A lp coherent phonons were also found to 
be positively chirped with large oscillation amplitudes of ~ 0.13 A [12]. Further, collapse 
and revival of chirped coherent phonon oscillations were observed in the reflectivity data 
when the carrier density was increased further bey ond a critical carrier density levels (3.5 x 
10 21 cm~ 3 at 10K and 5 x 10 21 cm -3 at 296K) [13]. This behavior was explained in terms of 
dynamics of a phonon wave packet in an anharmonic potential, where the packet periodically 
breaks up and revives to its original form, implying a nonclassical dynamics. A detailed first 
principle density functional calculations and optical double pump pulse experiments show 
that the anharmonic contribution to the phonon period is negligible under such high pump 
fluence regime in bismuth [14] and tellurium JfJ 10| and time dependence of carrier plasma 
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density alone is able to explain the observed softening of the phonons. 

Now, we turn to Bi2Te3-subject matter of the present study. Bi2Te3 is an important 
material both from t he p oint of view of thermoelectric flf] ] as well as exotic physics of 



topological insulators 



17, 



181 ] . There have been two earlier studies 19|, |20j of Bi 2 Te 3 using 
non-degenerate pump probe experiments at room temperature with 100 fs pulses where 
pump ffuence was kept below 1 mJ/cm 2 . In ref 19|, the two observed coherent phonon 
modes at 1.85 THz and 3.68 THz were assigned as A\ g and its second harmonic. In a 
later study [20], the modes seen at 1.85 THz and 4.02 THz were assigned to two allowed 
Ai g modes, in agreement with conventional Raman measurements [21]. The objective of 
the present work is to study coherent phonons in single crystal of Bi 2 Te 3 as a function of 
temperature (from 296 K to 3 K) at high photoexcited carrier densities, with a view (i) to 
understand the assignment of the high frequency coherent mode (first order vis a vis second 
order) and (ii) to study phonon chirping. The Gabor wavelet transform has been performed 
on time domain data to study the evolution of the coherent phonon modes. We note that 



earlier studies on Bi 2 Te 3 [19 
(with negligible chirping). 



20] were done only at room temperature and low pump fluences 



EXPERIMENTAL DETAILS 

A single crystal of Bi 2 Te 3 (6 x 6 x 0.5 mm 3 ) with a cleaved surface perpendicular to the 
trigonal axis mounted on a continuous helium flow cryostat was used in our experiments. 
Femtosecond pulses were derived from Ti:Sapphire amplifier (Spitfire, Spectra Physics Inc) 
producing ~ 50 fs pulses with the central photon energy of 1.57 eV at a repetition rate of 
1 KHz. The pump beam was modulated at 393 Hz with a chopper and the reflected probe 
intensity was recorded using a Si-PIN diode and a lock-in amplifier. The spot size (half 
width at 1/e maximum) of the pump and probe beams were kept at ~ 600 fim and 400 /zm, 
respectively at the overlap of the two beams on the sample. The pulse width was measured 
to be 65 fs (full width at half maximum) using a thin beta-barium borate (BBO) crystal 
at the sample point. The polarization of the pump beam was kept perpendicular to that 
of the probe beam to avoid scattered pump light reaching the detector. Both the pump 
and probe beams were kept close to normal incidence. The crystal surface is seen to get 
damaged at pump fluence beyond ~ 4.5 mJ/cm 2 and hence all our experiments were done 
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at pump fluences of 3.3 mJ/cm 2 and 1.3 mJ/cm 2 , whereas the probe fluence was kept at 0.4 
mJ/cm 2 . The time resolved reflectivity of the sample was recorded as a function of sample 
temperature varying from 296K to 3K. 

Bi 2 Te 3 is a narrow band semiconductor with an indirect band gap of 0.15 eV and it 
crystalizes in the R3m structure with the point group D| d . It is made up of close-packed 
atomic layers which are periodically arranged along the c-axis in five layers (Te«-Bi-Te( 2 )-Bi- 
Te 1 - 1 )) called as 'quintuples'. These layers are bonded by van der Waals force and the weakest 
link among the layers is Te^-Te^. Five atoms per hexagonal unit cell (a = 4.38 A, c = 
30.49 A) give totally twelve optical phonons out of which four are Raman active modes 
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, |22| represented as 2 A l9 (observed at 1.86 THz and 4.02 THz) + 2 E 5 (observed at 1.1 
THz and 3.09 THz). The linear absorption coefficient, a of Bi2Te3 at 1.57eV is ~ 4 x 10 5 



cm 1 



23fl and thus the penetration depth, £ (~ 1/a) is 25 nm. This corresponds to PCD, 
iVo (= 6.25 x 10 18 Fa{\ - R)/E p ) ~ 1.7 x 10 21 cm~ 3 at pump fluence F = 3.3 mJ/cm 2 
at the sample surface, which is about 1% of all the valence electrons. Here R=0.68 is the 
reflectivity coefficient of Bi 2 Te 3 , at photon energy E p of 1.57 eV. 

RESULTS AND DISCUSSION 




5 10 15 
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FIG. 1: Normalized time resolved reflectivity change (AR/R) of bismuth telluride as the function 
of the time delay between the pump and probe pulses at various temperatures (a) T = 296 K, (b) 
150 K , (c) 90 K, (d) 10 K, (e) 3 K (colour online). 

The coherent phonon mode's normal coordinate is written as [24], Q = 
b exp(—7fyt)cos(27ii/ t+(f)) where b, 7, v and (ft are the amplitude, damping constant, fre- 
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quency and the initial phase of the coherent phonon mode. To first order in Q, the normalized 
change in reflectivity (AR/R) of a probe beam due to generation of two A\ g coherent phonon 
modes (A^ and A^) in an absorptive material can be written as 



AR _sr 

i=l,2 



d(AR/R) 
dQ 



bi exp(-TTjit) cos(2i\Vit + (pi) (1) 



Here i — 1 and 2 corresponds to A lg and A^ respectively. The time resolved reflectivity 
data for Bi 2 Te 3 at a few temperatures using pump fluence of 3.3 mJ/cm 2 are shown in 
Fig. [TJ The signal contains both non-oscillatory and oscillatory components. The non- 
oscillatory background arising from carrier dynamics was removed by a digital band pass 
filter to extract the oscillatory part jljllizOj. The oscillatory part of the transient normalized 
differential reflectivity data was analyzed using Eq. (Tj[|). The fit was satisfactory at room 
temperature but was not good at lower temperatures. This necessitated the inclusion of 
chirp parameter f3 in Eq. p]) as fl3| 

— = J2 B i exp(-7T 74 t) cos(2irVit + fat 2 + 0,) (2) 

i=l,2 

where Bi = [ ^^q ] bi is the coherent phonon amplitude and is the chirp parameter. 
This fit was found to be excellent over the entire temperature range. For example, Fig. 

displays the digital band pass filtered time domain data at 3K recorded using pump 
fluence of 3.3 mJ/cm 2 along with the fit using Eq. ([T|) (dashed line) and (T5J) (solid line), 
and Fig. [2] (b) shows the corresponding fast Fourier transform (FFT) of the data and the 
fits. Here the FFT intensity of the second mode in the frequency range of 3 THz - 5 THz 
is appropriately scaled up to compare with the first mode. The eigen vectors corresponding 
to the Aig modes [2j| are also shown as the inset of Fig. E(b). It may be very difficult 
to see the chirping effect in time domain as the higher frequency mode is very short lived 
(1.2 ps) and therefore, frequency domain is a better choice in identifying the chirp in such 
cases. The FFT of the fitted function in time domain (Eq. (T5])) with chirp (/3 ^ 0) shown by 
continuous line) fits the FFT of the measured data (open circles) much better than without 
(3, i.e f$ — (short dashed line) [Fig. [21(b)]. The value of /3 2 = 0.38ps~ 2 is much higher 
than 2 x 10~ 3 ps -2 ). When the pump fluence is decreased to 1.3 mJ/cm 2 , (3% = 0.045 

ps -2 which clearly indicates that the chirping is mainly due to large number of photoexcited 
carriers. Since our focus is mainly on the high density photo carriers mediated coherent 
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phonons, the temperature dependence of the fit parameters at highest pump fluence of 3.3 
mJ/cm 2 are discussed from here on. 




FIG. 2: (a) Digital band pass filtered normalized time resolved differential reflectivity data at T = 
3K (open circles) along with their fits according to Eq. ([T]) (dashed line) and (T2]) (solid line), (b) 
The corresponding FFT of the time domain data and the fits. The data in the range of 3 THz -5 
THz has been scaled appropriately to compare the asymmetry seen at 3K. Here, the short dashed 
line is the FFT of the fit without chirp parameter. The eigen vector of the two A\ g modes are 
shown in the FFT panel (colour online). 

The parameters (filled diamond) obtained from the fit of the time domain data are plotted 
as a function of temperature in Fig. [3J Though the dependence of f3\ and on temperature 
are similar, the chirp parameter @2 for the high frequency phonon mode is almost 

(2) 

two orders higher compared to (3\ showing that electron-phonon interaction for A\ should 
be much higher than that for A^ (Figs. [3(a) and (b)). The increased phonon chirping 
seen at low temperatures and its inverse dependence on temperature /3 ~ T _1 is shown 
using continuous lines. An increase of with decreasing temperature can be qualitatively 
understood by considering the carrier diffusion across the penetration depth of £ in Bi 2 Te3. 
At very high carrier densities (~ 10 21 ), it has been shown that the carrier diffusion coefficient, 



D a has major contributions from carrier-carrier scattering mediated diffusion, D e h |25l. |26|. 
In such a case D a ~ T 5 / 2 and hence the diffusion time t^h = £ 2 /D a ~ T~ 5 / 2 . This inverse 
power law dependence suggests that at low temperatures, the diffusion time is longer and 
hence electron phonon interaction will result in large chirping. The temperature dependence 
of chirping as /3 ~ T _1 needs new theoretical inputs. 
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Next, we turn to the phonon amplitudes. It is seen that both the coherent phonon 
amplitudes, Bi and B2 increase as the temperature is lowered (Figs. [31(c) and (d)). It 
has been shown |27j | that the coherent phonon amplitude, behaves quite similar to the 



27] since the source of spontaneous 



temperature dependence of the Raman peak intensity 
Raman scattering and the driving force in the generation of coherent phonons through ISRS 
are the same 28] . The Raman cross-section increases with temperature as \n{y) + 1] (where 



n(z/) is the Bose-Einstein statistical factor) and hence the Raman peak intensity can be 
written as, I p ~ [2^/2) +1] ' wnere the temperature dependence in the denominator comes 
from the cubic anharmonic contributions to the linewidth. The fit (solid line) using this 
expression is shown along with the data for Bi and B2 in Fig. [3]^c) and (d). 




FIG. 3: Phonon chirp parameters [(a) and (b)], amplitude of oscillations [(c) and (d)], damping 
constant [(e) and (f)] and frequency [(g) and (h)] of both the coherent Ai 5 modes obtained from 
time domain fit to the data with high pump intensity (Filled diamonds) versus sample temperature. 
The continuous lines in (a)-(h) are the fits (see text) (colour online). 

The temperature dependence of coherent phonon damping factor and the frequency of 
the two modes are displayed in Fig. E3e),(g) and Fig. [3]^f),(h), respectively. In the case of 
A-ig , the frequency decreases from 1.9 THz to 1.84 THz, i.e., a decrease of 3% as the crystal 
is heated from 3K to 296K and the damping term increases by 150 % from ~ 0.04 THz 
at 3K to ~ 0.1 THz at 296K. The behavior of the second mode is rather interesting: the 
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damping constant becomes almost constant after 60K and the frequency decreases from 4.0 
THz to 3.6 THz (10% change) when the crystal temperature increases from 3K to 296K. 



The fit (thick line) shown in the figure is by using the well known functions 29] based on 
cubic anharmonicity where the phonon of frequency, v decays into two phonons of equal 
frequency: J p h(T) = 7o + C[l + 2n(t/ /2)] and v p h(T) = + A[l + 2n(z/ /2)] where u 
(frequency at T=0K), A, C and 70 (disorder induced damping) are the fitting parameters 
(A and C are the measures of third order cubic anharmonicity). The parameters obtained 
from fitting are, vq = 1.908 THz, A = -0.005 THz, 7o = 0.040 THz and C = 0.005 THz for 
A^ mode; and v = 4.007 THz, A = -0.060 THz, 7o = 0.263 THz and C = 0.017 THz for 

(2) 

A\ mode. Thus the disorder induced damping (70) and anharmonicity (A and C) are more 
for A^g compared to A^j . It can be seen from Fig. E^li) that the temperature dependence 
of V2 is non-monotonic. The decrease of at 3K may be due to the chirp in the frequency 
whose quantitative understanding is lacking. 

Finally, to examine whether the second mode at ~ 4 THz is the second harmonic of the 

n 

first Ai s mode at 1.8 THz as suggested in ref [19], we re-fitted the digitally band pass filtered 
time domain data with ^ expressed in a second order approximation in Q as 19 1 



AR 
~R~ 



d(AR/R) 
dQ 



d 2 (AR/R) 
dQ 2 



Q 2 (3) 



where Q is for the A^ phonon mode. Taking the expression for Q as in Eq. flTJ (since /3± is 
small), the fit according to Eq. (j3J) was unsatisfactory as the second mode could never be 
accommodated in the fit to the data as shown in Fig. HJ In Fig. HJ the time domain data 
(thin line) along with the new fit (thick line) using Eq. (j3J) and their FFT are shown in (a)- 
(b) for 3K and in (c)-(d) for 296K where open circles are the FFT of the data and thick lines 
are the FFT of the fit. It can be seen from these figures that the fit is not compatible with 
the data in both frequency and time domain. At 296K, though the fit seemingly agrees with 
the data in the frequency domain, it is not so in time domain. This confirms that the higher 
energy phonon mode is not a second order mode of A^. To further elucidate this point, 
we have analyzed the temperature dependence of the ratio of the integrated intensities of 
coherent phonons, (open circles in FigJU^e) where the y-axis is normalized to the ratio 
at 3K). Here, we consider two cases: (i) For V\ and v 2 with v 2 > v\ as the single phonon 
modes (SPM), the ratio of their integrated intensity, is ~ ^|^)+i shown by solid line in Fig. 
Hie), (ii) for z/ 2 as the second harmonic or two phonon mode (TPM) of V\ (i/ 2 = 2ui), we note 
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that the ratio of their integrated intensity is [n(v\) + l] 2 /[n(z/i) + 1] = n{i>\) + 1 shown by 
dashed line in Fig. H] (e). Thus, it is clear from Fig. H] (e) that SPM is compatible with the 
data corroborating our conclusion that the high frequency mode is a single phonon mode 
and not a second harmonic of A\„. 




FIG. 4: The time domain data (thin line) along with the fit (thick line) according to Eq. ([3]) and 
their FFT are shown in (a)-(b) for T=3K and (c)-(d) for T=296K. In (b) and (d), the open circles 
are the FFT of the data and thick line is the FFT of the fit. The ratio of integrated coherent 
phonon intensity (open circles), along with the fit using single phonon (SPM) (thick line) and 
two phonon model (TPM) (dashed line) are shown in (e). The estimated lattice displacements for 
the two coherent phonon modes (open star for Ui and open inverted triangle for U2) are given in 
(f) along with a linear fit (colour online). 



The lattice displacement of the coherent pho non modes at 1.9 THz and 4.02 THz can be 
estimated for absorbing materials using Jill. I30I 



Uf 



3.8 x irr 3 Bi F 

QUi\e\ 



' 2eaV 

D 



(4) 
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where Z7» is in Angstrom (A), F is the pump fluence in mJ/cm 2 , g is the density of the 
material in amu/A 3 and e is the dielectric constant (e\ + jsz), ^ P h is the energy of the 
phonon in eV and D = -^||| with E as the photon energy in eV. Now, for BioT^, D = 
iff ~ lO^eV- 1 at 1.57 eV and |g ~ 10 3 eV~ l with e = 2.75 + j 15 at 1.57 eV |23||. Thus, 
for Bi 2 Te 3 , the lattice displacement is 



38 F 

U^^B,—- (5) 

QVi\e\ 

The temperature dependence of Ui and U2 thus estimated is shown in Fig. H] (f), which 
essentially arises from temperature dependence of Bj and z/j. A linear fit (line) to U is shown 
in the figure. 

To capture the evolution of coherent phonons with time, we have performed continuous 
wavelet transform (CWT) similar to that used by Hase et a 
a quasiparticle in silicon. We have used the MATLAB code 
wavelet based on Gaussian function given as 33^ - 



^_ 3l| in observing the birth of 



32| modified for Gabor mother 



1 \ 1/2 

m/s) = tt" 1/4 ( — exp 
\psj 



t 2 t 
+]- 



2s 2 p 2 



(6) 



where s is the scaling factor (inverse of frequency) and p=n(2/ln(2)) 1 / 2 is a constant. Here, 
we describe the procedure to calculate the continuous wavelet transform. The wavelet trans- 
form of a given time signal, x(t) is given by, 

CWT(t,s) = — L / x(t)^*((t-r)/s)dt (7) 



A starting scale, s s t ar t corresponding to a frequency higher than the highest frequency of 
the signal (determined by FFT) is chosen and the starting wavelet \& s tart{t I ' s start) i n time 
domain here is a compressed wavelet. The cross correlation of ^ s tart{t) with x(t) is computed 
using Eq. with r = 0. The magnitude of cross correlation will depend on how closely the 
frequency components in x(t) and ^ s tart{t) match. This procedure is repeated by translating 
^start{t) in time domain (r) and this gives wavelet coefficients for a given s and a range of 
values of r. The whole procedure is again repeated for next higher s. Thus, one gets a range 
of wavelet coefficients of x(t) in time-scale plane which is then converted into time-frequency 
plane. The wavelet transform of a time domain signal gives three dimensional (3D) plot of 
wavelet coefficients vs frequency and time. CWT chronograms (contour of the 3D plot) of 
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FIG. 5: (a) CWT chronogram of high pump intensity data at T=3K. The color bars are given at the 
right side to symbolically represent the CWT intensity, (b) The maximum of CWT intensity at each 
time delay for both modes (thicker line is for mode and the thinner line is for the mode). 
The time delay between the two modes (open stars) are plotted as a function of sample temperature 
in the inset where the thick line is the guide to the eye. (c) The frequency corresponding to the 

(2) 

maxima of CWT intensity at each time delay for A\ g ' (thicker line is the linear fit to the data) and 
for (d) A^ mode, (see text) (color online). 

the time domain data at 3K and high pump excitation (3.3 mJ/cm 2 ), are shown in Fig|5](a). 
It can be seen from Fig [5] (a) that at 3K, the frequency of mode starts at 1.5 THz and 
reaches 1.9 THz in ~ 750 fs and this time is seen to be constant at all the temperatures. We 
could not resolve a similar build up in frequency for the A^ due to its short lifetime (0.9-1.6 
ps) compared to the first mode (4-10 ps). To get more insight into these chronograms, the 
maximum of CWT intensity at each time delay (time slice) was calculated and plotted at 
3K for both the modes (Fig. |5^b)). Here, thicker line is for A^ and the thinner line is for 
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AL . The time to reach the maximum of CWT is ~750 ± 20 fs fs for A\ mode and ~ 510 
± 20 fs for A\ g mode. The time delay of ~ 240 ± 20 fs between the maxima of the two 
modes A\ g and A\ g is found to dependent on the sample temperature as shown in the inset 
of Fig. |5] (b) . The reason for the time delay and its temperature dependence is yet to be 
understood. The phonon chirp for the second mode is demonstrated through the plot of 
frequency(open circles) corresponding to the maxima of CWT intensity at each time delay 
at 3K in Fig. [5] (c). Thicker line is the linear fit to the data with v 2 {THz) = 3.8 + j3 2 / (2n)t 
where t is the time delay in ps. This corresponds to 02 = 0.41 ps~ 2 , closely matching to 
the value of 0.38 ps -2 derived from the time domain fit. The oscillatory structures seen in 
the frequency is an artifact in the wavelet transform when the separation between the two 
frequencies is less than 4 THz. Similar analysis for A\ mode is given in Fig. [5] (d) where 
the build-up time is ~ 750 fs (consistent with Fig. [5] (b)) and the frequency is seen to be 
constant with delay time as exp ected since (3\ is very small. We note that the above analysis 
using AGU-Vallen-wavelet [34] resulted in similar results. 



CONCLUSIONS 



In conclusion, we have studied coherent Ai g phonons in Bi2Te3 as a function of both 
temperature and pump fluence. We have observed that the higher frequency coherent phonon 
mode at ~ 4.0 THz is and is not a second harmonic of A^. It acquires two orders 
of magnitude higher positive chirping at the lowest temperatures and high pump fluence as 
compared to the A^ phonon. The chirp in time domain is manifested as an asymmetry 
in the frequency domain. The wavelet transform of the time domain differential reflectivity 
helps to identify the chirp and the time for a phonon mode to build up to its maximum 
amplitude. We hope our experiments will motivate theoretical calculations to understand 
the exact temperature dependence of phonon chirping at highly excited carrier densities, 

(2) 

higher chirping for A lg and the different build up times for two coherent modes. 
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